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Abstract 

We consider the variant of stochastic homogenization theory introduced in [71 [8]. The equation under con- 
sideration is a standard linear elliptic equation in divergence form, where the highly oscillatory coefficient is the 
composition of a periodic matrix with a stochastic diffeomorphism. The homogenized limit of this problem has been 
identified in [Jj. 

We first establish, in the one-dimensional case, a convergence result (with an explicit rate) on the residual process, 
defined as the difference between the solution to the highly oscillatory problem and the solution to the homogenized 
problem. 

We next return to the multidimensional situation. As often in random homogenization, the homogenized matrix 
is defined from a so-called corrector function, which is the solution to a problem set on the entire space. We describe 
and prove the almost sure convergence of an approximation strategy based on truncated versions of the corrector 
problem. 

7— I 1 

>: 

^ ; 1 Introduction 

^-j- ■ Homogenization theory for linear second-order elliptic equations with highly oscillatory coefficients is a well developed 
• \ topic. In the periodic case, the homogenized problem is known, and convergence rates of the oscillatory solution 
I ■ (denoted u e ) towards the homogenized solution u* have been obtained. 

\ The situation is less clear in the random (say stationary ergodic) setting. The convergence of u e {-,uj) to some 
y—i • deterministic it* is a classical result. However, rates of convergence are much more difficult to obtain. A central 

\ difficulty in stochastic homogenization is that the corrector problem, that needs to be solved to next compute the 
• '-h ■ homogenized matrix, is set on the entire space (in contrast with the periodic case, where it is set on the periodic cell). 

[ This induces many theoretical and practical difficulties. 
^ In what follows, we are interested in the problem 



div 



(^,u?jVu e (x,Lu) = f(x)w.V, u e (-,Lu) = on dV, (1) 



where the random matrix A satisfies standard coercivity and boundedness properties (and some structure assumptions 
that we detail below), T> is an open bounded set of R d and / e L 2 (V). 

The analysis of the residual, that we define as the difference between the oscillatory solution u £ and the homog- 
enized solution u±, was first taken up in [5], and next complemented in [3!. Both studies consider the equation 
du e ~\ 

f(x) in the one-dimensional setting, where a (x,ui) is a random stationary process. The behav- 



d_ 

dx 



fx \ au 
\ e ' J dx 



ior, when e — > 0, of the residual u e {x,uj) — u*(x) turns out to depend on the asymptotic behavior of the correlation 
function of the conductivity coefficient rj{x) := Cov(a(0, •), a(x, •)). In [5], the case of small correlation lengths is 
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studied, which amounts to assuming that rj{x) ~ x ->oo % a with a > 1. The correlation function is thus integrable. In 

u £ (x, lu) — uJx) „ 
that case, when e — > 0, the random process = converges m distribution to a Gaussian random process. 

The case of long correlation lengths, namely when rj(x) ~x->-oo % a for some < a < 1, is studied in [3], where it 

u £ (x, cu) — Ui,{x) 

is shown that the random process converges in distribution to a Gaussian random process (defined 

using a fractional Brownian motion) . This result shows that the rate of convergence of u £ to can be as slow as e Q / 2 
for any a > 0, without any further assumptions on the stationary process a. 

Of course, in both works, the one-dimensional setting allows to get some analytical expression for the residual. In 
turn, the analysis of the asymptotic behavior of the residual performed in [3] relies on this analytical expression. In 

higher dimensions, the case of the equation — Au £ + q ^ — ,tuj u £ = f(x) has been studied in [2]. 

Our first aim here is to study a similar question for a variant of the classical stochastic homogcnization theory. We 
consider in the sequel the following problem, which has been introduced in [7j and further studied in [8]: 



div 



-4, 



per (j)' 1 ^u £ {x,uj) =f(x)mV, u £ {-,u) = on <9X>, (2) 



where <j) is almost surely a diffeomorphism from R d to R d with some stationary properties and A pcl is a Z d -periodic 

matrix, that satisfies the classical coercivity and boundedness properties (see precise assumptions in Section [2J] below). 

This model is appropriate to represent a periodic, ideal material, that is randomly deformed (think of fibers in a 

composite material that are placed at a random position, rather than on a perfect, periodic lattice). In 7 ], it is 

shown that the solution u £ (-,u) to the above problem converges as e goes to to it*, solution to some homogenized 

problem (see Section [5] below) . In the sequel, we aim at obtaining the rate of convergence of u £ to u*, in the one 

dimensional setting. We make below an assumption on the random diffeomorphism which implies that our setting is 

close to the one studied in [9 (rather than that studied in [3j). Under this assumption, we show that the random 

u £ (x, uj) — u+fx) „ , 

process = converges m distribution to a Gaussian random process that we completely characterize (see 

Section EH Theorem H). 

We next turn to a question of different a nature. As pointed out above, the homogenized matrix A* associated to © 
depends on the solution of the corrector problem, which is set on the entire space. Computing an approximation of A* 
is thus, in practice, a challenging question. A standard strategy is to consider the corrector problem on a large, but 
bounded domain Qn, supplemented with (say periodic) boundary conditions. An approximation of the exact corrector 
is thus computed, from which an approximate homogenized matrix A* n (uj) is inferred. As a by-product of working on 
a bounded domain, the approximate homogenized matrix is random. In the classical random homogenization setting 
(that is ((T|) where A is a stationary matrix), the convergence (and its rate) of A^(uj) to A* has been studied in [15] . 
using some previous approximation results [35J. It is shown there that ^^(a;) almost surely converges to A*, and 



that E 



\A* N -A*\ 2 



converges to as N a , for some a > which implicitly depends on the mixing properties of the 

random coefficient A of the equation (fT]). It is expected that, depending on the properties of that random coefficient, 
a can be arbitrary small. 

In this work, we consider the above variant ([2]) of the classical random homogenization setting. We describe 
a strategy (originally introduced in 12 ) to approximate A* which is based, as in the classical setting, on solving 
the corrector problems on bounded domains Qn- We prove here the convergence of this approach (see Section EOI 
Theorem g]) . 

Our article is articulated as follows. In Section [5] we present in details the variant of the classical random ho- 
mogenization introduced in pj [8]. We next present in Section |3] our two main results, first on the residual process in 
dimension one (see Section 13.11 and Theorem [5]) , second on a practical approximation of the homogenized matrix in 
dimension d > 2 (see Section I3T21 and Theorem 2]). The subsequent two sections are devoted to the proof of Theorem^ 
The actual proof is performed in Section 3] and needs some technical results which are proved in Section [5] Our final 
section, Section [5] collects the proof of Theorem 0] 
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2 A variant of the classical random homogenization 



To begin with, we introduce the basic setting of stochastic homogenization we will employ. We refer to [T3] for a general, 
numerically oriented presentation, and to [TJJ [TJ] for classical textbooks. We also refer to [7J |S] for a presentation of 

our particular setting. Throughout this article, (f2, J 7 , P) is a probability space and we denote by E(A) = / X (uj)dP(u)) 

the expectation value of any random variable X E L 1 ^, dF). For any fixed d E N* (the ambient physical dimension), 
we assume that the group (Z d ,+) acts on il. We denote by (Tk)kez d this action, and assume that it preserves the 
measure P, that is, for all fc G 7L d and all B E T , FfaB) = P(B). We assume that the action r is ergodic, that is, 
if B E T is such that t^B = B for any k E Z d , then ¥(B) = or 1. In addition, we define the following notion of 
(discrete) stationarity (see [TJ El) : any F E L\ oc (R d ,L 1 (fi)) is said to be stationary if 

Vfc E 7L d , F{x + k,uj) = F(x,TkUj) almost everywhere and almost surely. (3) 

In this setting, the ergodic theorem [TBI HH1 HI] can be stated as follows: Let F E L°° (M d ,i 1 (il)) be a stationary 
random variable in the above sense. For k — (ki, k^ 1 ■ ■ ■ , kj) E 1, d , we set \k\oo = sup \ki\. Then 

Ki<d 



lW V F(x,T k oj) — > E(F(x,-j) in L°°{R d ), almost surely. 

(2N + 1)" * — ' N-¥oo 

|fc|oo:<iV 

This implies that (denoting by Q = (0, l) d the unit cube in Mr) 

F^,w)-^W.\^j F(x,-)dxj inL°°(R d ), almost surely. 



2.1 Mathematical setting and homogenization result 

As pointed out in the introduction, we consider in this article the following problem, which has been introduced in [7] 
and further studied in [5]: 

A pcr (V 1 (f > w )) V!i £ (a),w)] = f(x) in V, u £ (;uj) = on dV, (4) 

where V is a bounded open set of M. d , f E L 2 (V), <fi is almost surely a diffeomorphism from M. d to R d , and A pcr is a 
Z d -periodic matrix, that satisfies the classical coercivity and boundedness properties: there exists a + > a_ > such 
that 

V£ e R d , a-|£| 2 < A peT (x)£ ■ £ almost everywhere on R d , and a + = \\ A pcr \\ Loo (^d^ < oo. (5) 
In addition, we assume that the map </>(•, uj) satisfies 

Esslnf (det(V0(x,w))) = v > 0, (6) 

EssSup |V0(x,w)| = M < +oo, (7) 

ojEQ, x£R d 

V</> is stationary in the sense of (8) 

Assumptions ([5]) and ([7]) mean that is a well-behaved diffeomorphism, uniformly in ui. Note that A per o (j)^ 1 is in 
general not stationary. The above setting is thus not a particular case of the classical stationary setting. 

In [7], it is shown that, under the above conditions, u e (-,u;) converges to u* almost surely (strongly in L 2 (T>) and 
weakly in iJ 1 (2?)) when e goes to 0, where u+ is the solution to the homogenized problem 

— div [A*Vu*(a;)] = f(x) in T>, = on <9X>. (9) 

In ([5]), the homogenized matrix coefficient A* is equal to 

VI < i,j < d, A*j = det (e Qf V<j>(y, -)dy^ E (j ) efA pm . {y, •)) (e, + Vw ej (y, •)) dyj , (10) 
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where Q = (0, l) d and where, for all p £ M. d , w p solves the following corrector problem: 



' -div [A peI (^(y.w)) ( P + Vw p (y,cj))] = in R d , 
w p (y,oj) = w p ((f)~ (y,u>),u), Vu> p is stationary in the sense of ©, 



(11) 



E| / Vw p (y,-)dy) = 0. 
'HQ,-) J 



2.2 The one-dimensional case 

Our first main result, presented in Section 13.11 is a convergence result in the one-dimensional case. In that setting, it 
is possible to write some explicit formulas. Choosing T> = (0, 1), the problems (jlj) and © respectively read 



d 
dx 



'[>■ i- \ 1 ' \ 71 UJ ) ) ~~fa( x,UJ ^ 



= f(x) in (0, 1), U £ (0, W ) = 0, u £ (l, w ) = 0, 



and 



The corrector problem (fTTj) . that reads 



aw 



-!(!,,«)) (l + —( V ,a;) 



in R. 



dw 



w(y,w) = w((f> 1 (y 7 u>),uj), —— is stationary in the sense of ([3]), 

dy 



e(7 ^(l/,W)=0, 



can be analytically solved. Its solution w satisfies 



dw . . a* 

1 + -ryy^) = — tz-u — i 



where the homogenized coefficient a* is given by 

(a*)" 1 = 



E 



1 



dy 



(12) 
(13) 



(14) 



(15) 



(16) 



As pointed out in [TJ, we observe on (IT51) that, in the one-dimensional case, the gradient of the corrector has the same 
structure as the highly oscillatory coefficient in (fT2")l : it is equal to a periodic function composed with (j)^ 1 . This is not 
the case in dimensions d > 2, as shown in [Jj. 

3 Main results 

In this article, we show the following two main results, Theorems [5] and [U 



3.1 Residual process in dimension one 

Our first aim is to characterize how the residual process u £ (x,u>) — u±(x) converges to zero, where u £ solves (fT2|) and 
u* solves (|13|) . To this aim, we make the following assumptions. Let us introduce the 1-periodic function 



tp(x) 



1 1 



a pcr (x) a* 



and the random variables 



rk+l 

Y k (u) = / TP(t)4>'(t,Lj)dt. 



(17) 



(18) 
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As ip is periodic and <f>' is stationary, the random variables Yj~ are identically distributed. Due to (jTHJ) , we have 

M(Y q )=m(J ^(t)0'(t, =0. 

We furthermore assume that the random variables Yfc are independent, and hence that 

the variables Y k are i.i.d. (19) 

Likewise, we consider the random variables 

D k (uj) = / tf(t,u)dt, (20) 



which are identically distributed, and make the assumption that 

the variables D k are i.i.d. (21) 
Remark 1. Suppose that the derivative of the random diffeomorphism 4> reads 

<l>'(y>v) = i + ^2x k (uj) Gper(y) i[fe,fe+i)(y), 

feGZ 

where X)-(ui) are independent and identically distributed random variables and G pcr is a l-periodic bounded function, 
such that, for some < m < 1, 

|-^o(^)| < fn almost surely and || G pe r||i°°(K) < fn. 

Then, the conditions ((BJ), © an d © are satisfied with v = 1 — m 2 > and M = 1 + m 2 . By construction, the 
assumptions (|19[) and (|21[) are also fullfilled. 



The first main result of this article is the following theorem, the proof of which is postponed until Section 

Theorem 2. Assume that a pel and <f> satisfy ([5]), (|6]), ([7]) and (|8|). Assume furthermore the independence condi- 
tions (|19l) and (|21[) . We consider u e solution to (I12[) and it* solution to (|13l) . Then the residual process converges in 
distribution to a Gaussian process, 

u e (x,u) - u*(x) c „ , s 
7= > G (x,u;), 

where 

Go(z,^) = ^E5L / K (x,i)dW t , (22) 



where Wt denotes the classical Brownian motion and Ko(x,t) is given by 

K (x,t) = (l M (i)-i) F(s)ds-F(tfj with F(t) = J f(s)ds. (23) 

Remark 3. it might be possible to weaken assumptions (|19p and (|2ip . and to on/?/ assume that the identically distributed 
variables are swc/i i/iat |Cov(Y"o? ^fe)| < +oo, and likewise for D k - We have however not pursued in that direction. 

fcez 

3.2 Approximation of the homogenized matrix 

In this section, we return to the multidimensional setting. To compute the homogenized matrix A* defined by (|10[) . 
we first need to solve the corrector problem (jlip . which is set on the entire space. In practice, approximations are 
therefore in order. 

In the sequel, we describe a strategy introduced in |12j . and that mimicks the approach proposed in [10 to approx- 
imate standard corrector problems in classical random homogenization. In this article, we analyze this approach and 
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prove its convergence (see Theorem 2] below) . This is our second main result. We refer to [TJ Section 3.2] for some 
illustrative numerical tests. 

Convention: Following [7J Lemme 2.1], we adopt the convention that [V^- = 7— - for any 1 < i,j < d. Hence, 

for any scalar-valued function tp, the gradient of tp — tp o (p is given by Vtp(z) = {V<j)(z)) T S7tp((p(z)). This convention 
implies that [Vc^" 1 )] V(0~ 1 ) = Id. 

Presentation of the approximation The weak formulation of the corrector problem ([TTj) reads as follows (see [7]): 
for all tp stationary in the sense of , we have 



E / (VtP(y,u)) T A peI (<p-\y, oj)) (p + Vw p {y,u J )) dy = 0, 
\J<KQ,-) j 

where tp = tp o cp . The above expression can be rewritten, after a change of variables, as 



E 



det(V0) (vVf (V0)-M pcr (p + (y<t>Y 



= 0. 



Since tp, V</>, A per and Vw p are stationary in the sense of ([31), the ergodic theorem yields 

lim / det(V0) (V^) T (V^)- 1 ^. (p + (V<p)- T Vw p ) = a.s. 

N->oo \Q N \ J Qn V / 

where Qn = NQ. For a fixed N, we now define the approximate corrector w p as the Qjv-periodic function satisfying: 

for all tp Qw-periodic, f det(V^) (v^) T (V0)" 1 A pcl (p + (V^^Vu^) = 0. (24) 
JQn \ J 

Note that w p is uniquely defined up to an additive constant. 

In turn, recall that A* is defined by (TIT))) . After a change of variables, we infer from that equation that, for any 
1 < z, j < c?, we have 

A*j = det (e Qf V0(y, .)dy) ) E Qf det(V<^(y, •)) ejA pei (y) (e 3 + (V0)^ T Vw £j (y, •)) 
The ergodic theorem yields 

A% = lim (det (-^- [ V<K;u)) r^- [ det(V0)ef A pCT (e. + (V0)- T V«y 

It is thus natural to approximate A* by the matrix A^(uj) defined by 

1 ' 



a.s. 



A* N ^) = det [ -^-J V0(-,«)) Z?.\U). (2--,) 



where the matrix B^f(u) is defined by, for any 1 < i,j < d, 

[B*^ («) = 17^7 / det (V0) ef A pcr ( ej + (V0)- T V^) 
\Qn\ Jq n v 

1 



|Qjv| J<P(Q n ,uj) 



A pci (<t>- l (v,u)) (e 3 - + Vu^(»,w)) d tf , (26) 



A?7„, , ,\ _ ~N(i-ll 



where, for any p G Mr, w p is defined by (|24l) and where 



w p (y,w) = w p ((j) {y,u),u). 

Note that, as is standard in stochastic homogenization, the approximation A* n (lu) is a random matrix, even though the 
exact homogenized matrix A* is deterministic. This is a by-product of working on the truncated domain Qn rather 
than R d . 
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Convergence of the approach We prove in Section [5] below the following convergence result: 

Theorem 4. Let (j) be a diffeomorphism satisfying ([6]), ([7]) and ([8]) ; and A pcr be a periodic matrix that satisfies the 
ellipticity condition ([5]). Then the random matrix A* n (lo) defined by (|25[) converges almost surely to the deterministic 
homogenized matrix A* defined by (jTUJ) when N — > oo. 



4 Asymptotic behavior of the residual 

The aim of this Section and of the next one is to prove our first main result, Theorem [2] Using the one dimensional 
setting, we first establish a "representation" formula for the residual (see Section |4~T1 Theorem [5]). Using this formula, 
we are next in position to study the asymptotic behavior of the residual when e — > (see Section |4~2]) . Section[5]collects 
the proofs of some technical results used in Sections 14.11 and 14.21 



4.1 Representation formulas 

The following technical result will be very useful in the sequel. Its proof is postponed until Section [54 



Lemma 5. Assume that a por and (j) satisfy ([5]), ©, ([7]) and ([5]). Assume furthermore the independence conditions (|19[) 
and (|2ip . For any < a < f3 < 1, and any A G L°°(a,/3) with A' G L 2 (a,j3), define the random variable 



1 



Am $ 



t 



w\\dt 



(27) 



where the function ip is defined by (I17p . For any p G N* , there exists a deterministic constant C p independent of A, e, 
a and f3, such that 



V £ >0, E[Z £ (a,/3,.) 2p ] <C p [((3-aT+e^/ 2 ] \\\A\\%, {a>fi) + (f3 - af \\A'\\% {oi 



The above result heuristically implies that the quantity J A(t)ip ^0 1 ^-,lj^J dt is of the order of \fe. 

We will show below a convergence result for the random variables Z £ (a, f3, u) (see Lemma 1 1 01 below) . The bound- 
edness result stated in the above lemma is however sufficient for now. Using it, we indeed prove the following theorem, 
which is a key ingredient to prove Theorem [2] 

Theorem 6. Assume that a per and cf> satisfy ([7]) and ([5]). Assume furthermore the independence condi- 

tions p9p and (|21j) . Let u e be the solution to (|12j) and it* be the solution to (|13D . Then 



u e (x,uj) — u±(x) — J Kq(x, t) t/j ( (j) 1 1 dt + r e (x, w), 



(28) 



where Kq is defined by (|23p . ij) is defined by (|17[) . and there exists a deterministic constant C independent of e such 
that, for any e > 0, 



sup E|r e (a;, -)| < Ce and E 

xG[0,l] 



'e\\L 2 {0,l) 



< Ce 2 



In addition, for any p G N* ; there exists a deterministic constant C' p independent of e such that 

~\r e (x, •) - r s (y, -)\ 2p ] < C p e 2 » J {x - y)*P + e*-V* 



VeG(0,l), V(x, y) G (0, iy 



E 



and 



VeG(0,l), V(x,y)€ (0,1) 2 , E [|r e (a;, •) — r e (y, -)| 2p ] < C p e p (x — y) 



2p 



(29) 



(30) 



(31) 



In view of Lemma [SJ the first term of the right-hand side of (|2"5)l is of the order of y/e. The term r E , which is of the 
order of e in view of (121?)) . is hence a higher-order term. The bounds ([3T)|) and (l3T|) will be useful below to show that 
some random process is tight (see Section I4T121 Theorem [5]). 

Using the same arguments, we show the following result: 
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Theorem 7. Assume that a pel and <p satisfy l[5]), (0 and ([5]). ^4ssiime furthermore the independence condi- 
tions (|19l) and (|2ip. Letu 6 be the solution to (|12p . u* &e £/ie solution to (|13|) , andw be the corrector, which solves (|14[) . 
XTien 



d 



Mia;, a; 



x = a 



-i (r 1 (~ ,«)) ^(^ (V 1 (^)) * 



where if) is defined by (|17l) . ifi is given by 

Ki(t) = a* (F(t) - ^ F(s)ds^ iwtfi F(i) = jf /(a) ds, 
and there exists a deterministic constant C independent of e such that, for all e > 0, 



E 



sup \r E (x,- 
xe[o,i] 



< Ce and E 



sup \r E (x, ■ 

xG[0,l] 



< Ce 2 . 



(33) 



(34) 



Again, in view of Lemma [5J the two first terms of the right-hand side of (|52"j) are of the order of y/e. The term r e , 
which is of the order of e, is hence a higher-order term. 

Remark 8. It is easy to deduce from (|32[) . using Lemma\5\ and (|34[) . i/iai i/iere exists a deterministic constant C 
independent of e such that 



E 



^u £ (x,uj) - u*(x) - ew ~^^) 



L 2 (0,1) 



< Ce. 



Likewise, we deduce from (|2"8")) , irnng Lemma [5] and (|29p . that 



E 



||u £ (x,a;) -u*(a)||i 3(0)1) 



< Ce. 



Using the expression (|54p below, we infer from (|35p and (|36p £/ia£ 



E 



r (x, w) - it* (at) - ew -^-( x ) 



ff!(0,l). 



< Ce. 



(35) 



(36) 



(37) 



We recover (in the one- dimensional situation) a classical result of homogenization: the corrector w allows to obtain a 
convergence result in the H 1 strong norm. We refer to Theorem 3] for a corresponding result in classical random 
homogenization (in the multidimensional setting). 

The proof of Theorems [6] and [7] are direct consequences of Lemma [5] and of the analytical expression of u £ and it*. 

Proof of Theorem^ Introduce F(x) = / f(t)dt. The solution to (TT2")) reads 

Jo 

i ,. r F(t) 



u £ (x, u>) — c € {uj) 



where 



o a pcr (0 Jo a P cr(0 1 



dt, 



(38) 



(39) 



o a P cr (</> 1 



dt 
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Likewise, the solution of the homogenized problem ([§]) is 



u*(a;) = c / — - 

a* J a* 



dt.. 



where a* is given by (fT6f and 



F(t)dt. 



Step 1: Representation formula 

We compute the residual process using (HOI) and ([5g|l: 



u E (a;, w) — = c £ (oj) 



-tit — c 



(c £ (c) - c*) j%U~ l (1> U )) dt + ( c » - 



where ?/» is defined by (flTl) . We also infer from ((39|) that 

/ , v -1 

1 



c e (w) -c* 



a P cr (0 1 
1 



(F(t)-c*) 



1 



df 



/o a per (^i- 1 (f,w)) 

y 1 1 

where we have used that, in view of (141 p. we have / f-FYt) — c*) — di = 0. Observe now that 



/O "per 1 (|,^)) 

We then deduce from (gSJ) that 

c £ (w)-c* = a* / (F(t)-c*)^ 



where 



a pcr f 



^U" 1 (-,W 



(F(i)-c*)^ U" 1 (-,w 



Collecting g2]) and (JUJ, we write 







(J 


-)) 






-1 













df 



r s (x 7 Lo)+J^ (c* - F(t))iJ U- 1 (^,cjjj dt+ (F(t)-c*)V^- x f 



r e (x,u)+ I K { X) t)i> [ 4>~ L [ -,u) )dt 



with 
and 



K (x,t) = (l [M (t)-a;) (c*-F(t)) 



(46) 



In view of (|4Tj) , we recover the expression f|23[) of if o • We thus have written the residual in the form ([28 

Step 2: Proof of the bound flU) 

We first bound p e (w). We infer from (|45|) that 



|p e (w)| < a + a* 



dt 



(F(t) -c*)i/) f 1 hit 



(47) 



Using the Cauchy Schwartz inequality, we deduce that 
E(|p £ |) < ea+a* 



\ 



E 



1 



1>[4> 



dt 



\ 



E 



(f(t)-c^ f 1 



Using Lemma [5] with p = 1, a = 0, j3 = 1, A(t) = 1 and A(t) = F(t) — c*, we obtain that there exists a constant C 
independent of e such that 

H\p e \)<Ce. (48) 



We also deduce from (1471) that, for any pgN*, 
E(V| 2p ) < (a+a*f p e 2p ^ 



E 



-== / V ( 0" 1 



4p\ 



\ 



E 



4,A 



Using again Lemma [5j we obtain that there exists a constant C p independent of e such that 

E(|p £ | 2p ) < C p e 2p . 



(49) 



Using the obtained bounds on p £ , we now estimate r e . We infer from (|44p . using (|4"9")l and Lemma [51 that, for any 



E 



(V-c^) < (a*) 2p e p 



C„E 



1 



(f(i)-c^U 



rff 



2p> 



C p E(V| 2p ) 



< (a*) 2p C p e p + C p e 



2p 



< C p e p 



(50) 



for a constant C p independent of e. In view of (14"6"1) . we thus obtain, using (|4"5)l and (|50p . that 



|r £ (z,.)|) < -E(|p £ |)+E 



(c £ - c*) / V 



< C e + ^E(|(c £ -c*)| 2 ) 



\ 



E 



— / ^ ( 



for a constant C independent from e and x £ (0, 1). This concludes the proof of the first assertion in (129")) . 
Similarly, we have 

2 



(r £ (x, u)f < tJ)IP» 2 + 2(c £ H - c*) 2 jf ( 
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thus 



E 



r ellL 2 (0,l) 



; Ce 2 + 2 



- [(c £ - c*) 4 ] E 



v^- 1 ) )<// 



Using (1501) and Lemma [S] with p = 2, we deduce that 



< Ce 2 + 2 / v^dx < Ce 2 
Jo 



l r e|lL 2 (0,l) 

for a constant C independent from e. This concludes the proof of the second assertion in (1291) 



Step 3: Proof of the bounds and (31]) 

We first prove (|30|) . In view of (|46|. we have 



r e (x,u)-r e (y,u)) = ^^p £ (uj) + (c £ {lj) - c*) / i/> 



thus 



and 



\r e (x,w) - r e {y,u)\ 2p < C p ( V (p») 2p + C p (c» - c*) 2p 



j>\ 0- 1 ( -,u ) )dt 



2p 



E 



|r £ (x,.)-r £ (y,-)| 2p < C p ( ^ ) E [(p e ) 2p ] + C p \/E[(c 6 — c*) 4p 



2)J 



E 



4/j' 



(» - xf P + y/ \x - yfP + £P-V2 



< C p £ 2p 

Since |y — x| < 1, we have (y — x) 2p < \y — x\ p < \J (x — y) 2p + sv^ 1 / 2 , and thus 

E 



\r e (x, •) - r s (y, -)\ 2p \ < C p e 2p ^(x - y) 2 P + 
This concludes the proof of (|3T)1) . We finally prove (T5TT) . We infer from ([5T|) that 



|r e (x,w) - r e (y,cj)| 2p < C p 



y-x 



2p 



(p») 2p + C p (c £ (u) - c*) 2p (x - y) 2p U\\ 



2p 

L°° (R) ' 



hence, using (gSJ and (1501) , 



E 



|r 8 (x, •) - r e (y, • )!* < C p (x - yf p [e 2p + e p ] 



This concludes the proof of (|3Tj) and thus that of Theorem [6l 



Proof of Theorem^ Recall that the solution to the corrector problem (|14p satisfies (fTS)) . We thus have, using 
and ([iDjl. 



d I fx \ du+ \ du" du* . . / , /a; \\ ortt* , fx 

-{u (x,u)-u*(x)-ew {_-,«)— (x)) = — {x ,u>)-—(x)(l + W (~,u>))-s^(x)w(- 



du 



dUi, 



d 2 u+ 



c £ (uj)-c* 



d 2 u± 



(x)w (^,^J 
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Using (|4"4"1) . we deduce that 

f™-o*(*-'(^))*-«£w.(H ^ 



di- 



aper (f 7 " V V £ ' // da;2 ' Opcr (</> -1 (f ,w)) 



with ifi defined by (15"3")) and 
Observe now that, in view of (|15l) and (IT7)) , we have 



r e M = ^ (53) 







(?») 






./() 



w(y,w) = a* / ^ (V^ 1 
Jo 

where we have chosen the integration constant in w such that w(0,u)) = almost surely. Thus 

) =a* jy ^(4>- l {t, u:))dt=^ j^L- 1 (l>")) dt - ( 54 ) 

Collecting this equation with ([52]) yields (J32J). The bound (JMJ) follows from ([531), (03) and (gSJ). This concludes the 
proof of Theorem [7] □ 

4.2 Proof of Theorem |2] 

T u e (ir, w) — uJx) ...... _ 

In this section, we prove that the random process = converges m distribution to a Gaussian random 

process that we characterize. Using (|2"5)l . we see that 

u e (x,uj) - u*(x) 

■= = G e (x,lj) + R e (x,cj), (55) 

where 



In view of (f2T?l) . we have 



G £ (x,w) = -j= / JT (a:,t)^ ( 0" 1 ( -,w ) ) A, (56) 
R e (x,u) = -pr e (i,w). (57) 



sup E|i? £ (x,-)| < C^e 
ce[o,i] 



for a constant C independent of e. As a consequence, 

Vx € (0, 1), i? e (a;, •) converges to in probability. (58) 

We are thus left with studying the behaviour of G e (x, lS) as e — s- 0. 

To prove that the random process G s (x, ui) converges in distribution, we will use the following result: 

Theorem 9 ([BJ, page 54). Suppose that (G £ ) £g ^ ^ and Gq are random processes with values in the space of continuous 
functions C (0, 1) with G e {0,u>) — Gq(0,cv) — almost surely. Assume that 

(i) for any k G N* and any < X\ < ■ ■ ■ < Xk < 1, the random variable {G e (x\,bS), . . . , G e (xk, w)) G converges 
in distribution to the random variable (Gq(xi,u>), . . . , Go(i£fc,u;)) as e — > 0. 
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(ii) {Ge) e( z(Qi) is a tight sequence of random processes in C (0, 1). A sufficient condition for the tightness of 
(G £ ) £g ^ i% is the Kolmogorov criterion: there exist 5 > 0, f3 > and C > such that 

We 6(0,1), V(x,y) G (0,1) 2 , E [\G e (x, •) - G £ (y, -)f] < C\x - y\ 1+s . (59) 

Then the process G e converges in distribution to the process Go as e goes to 0. 

For any x G (0, 1), the random variable G e (x,oj) is of the form of the random variable Z e (a, f3,uj) defined in (I27[) . 
with a = 0, f3 = 1 and A(t) = Ko(x,t). In Lemma[5j we have shown that the random variable Z £ (a, f3,uj) is bounded 
in the L 2p norm. We now show that this random variable converges in law to a Gaussian random variable. This will 
be a key ingredient to prove the first condition of Theorem [5J 

Lemma 10. Assume thata pcv and<p satisfy ([5]), © and ©. Assume furthermore the independence conditions (|19l) 
and (|21[) . For any < a < /3 < 1, consider a function A, piecewise continuous over {a, ft), with a finite number of 
discontinuities located at points {tk}i<k<m, and such that A' G L 1 ^, tfc + i) for any 1 < k < m — 1. Consider the 
random variable 

Z £ (a,(3,uj) = ±=j A(t)if> (V 1 d< ( 6 °) 

where the function ij) is defined by (JT7J) . Then Z s (a, converges in distribution to a Gaussian random variable 

Zo(a, (3, us), of mean zero and variance a(a, /3) = j °^ H^ll^fo, p\> withYar(Yo) = E 

E(/ ft) 



We write 



\MJqV) Ja 

where Wt denote the classical Brownian motion. 



Z (a,/3,u>) = ^ Yar(Xo) f A(t)dW t , (61) 



The proof of Lemma [TO] is postponed until Section [ 
To prove the second condition of Theorem [SI we will show that G e (x,uj) satisfies (|59l) . Observe that 



1 r , ( ,_i It \\ ,. £ Z" 1 ft 



with ,4(i) = / -F(s) ds — F(t), where F(t) = / f (s) ds. To prove that Gr(x. ui) satisfies ([5^|) . we will use the following 

Jo Jo 
result, the proof of which is postponed until Section 15.31 

Lemma 11. Assume thata per and(f> satisfy (O, © and ([5]). Assume furthermore the independence conditions Q19p 
and (|2"Tj) . Consider two functions A\ and Ai with Aj G L°°(0, 1) and A'j G L 2 (0, 1), j = 1,2. For any x G (0,1), 
consider the random variable 

H e (x,oj) = ^=£ Ai(t)iP L- 1 (^,u^dt+^=j\ 2 (t)^^f > - 1 (~,"j\dt (63) 

where the function ?p is defined by (JTTJ) . 

For any p G N* ; there exists a deterministic constant C p independent of e, x and y such that 

V £ G(0,1), V(x,y)e (0,1) 2 , E[\H £ (x,-)-H £ (y,-)\ 2p ] < C p (\x - y\? + e^ 1 " 2 ) . (64) 

In addition, there exists a deterministic constant C independent of e, x and y such that, for any x and y with \x — y\ < e, 



\H s (x,lj) -H e {y,u)\ < Cy/\x-y\ a.s. (65) 
We are now in position to prove Theorem [21 
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Proof of Theorem^ We have seen (see (|55p ) that 

u £ (x, to) — u ± {x) 



G e (x,uj) + R s (x,u>), 



(66) 



where G e (x,uj) and R e (x,uj) are defined by (|56")) and (|57|) . respectively. 
Let us study the process G £ (:e,cl>), which reads, we recall, 



G e {x,w) 



1 



K (x,t)fl> [cj)- L \-,u;) ) dt. 



As Ko(0,t) = for any t, we have that G e (0,w) = for any e, almost surely. We first show that this process 
satisfies the first condition of Theorem[9] For each set of points < x\ < ■ ■ ■ < Xk < 1 and each X = . . . , £ M. k , 
we consider the random variable 

k 
3 = 1 



Observing that 



we can write z F as 



where 



G e (xj,u) 



K {x ,t)^ U -1 - w eft 



A x (t)=J2^Ko(xj,t) = ( f F{ s )ds - F(t)) ^j(l [0<Xj] (t) 

,=! Vo /, = 1 



with -F(i) = / f(s)ds. By assumption, / € L 2 (0,1), thus _4;t is piecewise continuous with a finite number of 
Jo 

discontinuities located at {xj}i<j<k- In addition, we see that, over each (:Ei,2i+i), 1 < i < k — 1, 

E^-E&pW 

is in L 2 (xi,Xi+i) C L 1 (xi,Xi+i). Thus, using Lemma ITUI we obtain that z £ (lu) converges in law to 



where Go is defined by (PZ!?1) . This implies that 



lim ] 

e->0 



cxp ( i^Cj'G ! e (^ ) •) 



lim E(exp(iz e )) = E(cxp(izo)) = E 



exp [ i^2^jG (xj,- 

3 = 1 



Hence, for any fc 6 N* and any < x\ < ■ ■ ■ < xu < 1, 

(G e (xi, lu), . . . , G £ (xfc, cj)) converges in distribution to (Go(xi, w), . . . , Go(xk, u)) as e — » 0. 
Collecting (JHSJ), (JSTJ) and (|58"]i. we obtain that 

u e (x, w) — Ui,{x) 



the residual process 



satisfies Condition (i) of Theorem ® 



with the limit process Gq(x,uj) defined by (|22|) . 



(67) 



(68) 
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We now prove the Kolmogorov criterion, first on the random process G e (x,ui), next on the process =— -. This 

will show Condition (ii) of Theorem [§] Following (1521 , we write 

<?.(*,«) = ±[ Am (V 1 * -fj^m (V 1 (^)) * 

with .A(t) = / F(s) ds — F(t), where F(£) = / f(s)ds. The assumptions of Lemma QT] are satisfied, thus, for any 

Jo Jo 
peN*, there exists C p such that 

Ve e (0, 1), V(ar, y) € (0, l) 2 , E [|G e (a:, ■) - G e (y, -)| 2p ] < C P (|* - y\ p + eW 2 ) . (69) 
This directly implies that 

when \x-y\>e, E[\G e (x, ■) - G £ (y, -)\ 2p ] < C p \x ~ y^/ 2 . (70) 
When \x — y\ < e, using (|65|) . we see that there exists a deterministic constant C independent of e, x and y such that, 



E \G £ (x, •) — G s (y, -)\ 2p < C\x - y\ p < C\x - y\^-V/ 2 when \x-y\<e. 
Collecting ([70)1 and (|71D . we obtain that 

Vee(0,l), V(x,y)e(0,l) 2 , E [|G £ (x, •) - G £ (y, -)\ 2p ] < C\x - y\^' 2 . 
We now turn to the process R £ (x,lu). In view of ([57} and (|30p . there exists C p such that 

Ve e (0, 1), V(x, y) e (0, l) 2 , E [\R e (x, •) - i? £ (y, 0| 2p ] < C p e p y/ (x - y) 2 ? + e^l 2 . 
Hence, we deduce that 



E 



\R £ (x,-)-R £ (y,-)\ 2p ] < C p e p \x - y]^ 1 ^ 4 when \x - y\ > e. 



When a; — y\ < e, using (|3Tj) . we see that 



E 



\R £ (x, ■) - R £ (y, -)\ 2p \ < C p (x - y) 2p < C p e p \x - y\ p < C p e p \x - y^ 1 )/ 4 when \x-y\< 
Collecting ([71)1 and (|75p . we obtain that 

Ve e (0, 1), V(s, y) € (0, l) 2 , E [\R £ (x, •) - R e (y, -)\ 2p ] < Ce p \x - y^ 2 ^' 4 . 
We next write, using (|66|) . 



(71) 
(72) 

(73) 
(74) 
(75) 
(76) 



u e {x,us) - u e (y,u:) - u*(y) 



2p 



< C p \G e (x,u) - G £ (y, w)| 2p + C p |i? e (x, w) - Rs(y,u)r ■ (77) 



2/> 



Collecting (j?2")) and (|76p . we obtain that 
Ve 6(0,1), V0z,y)e (0,1) 2 , 1 
We thus obtain that 



u E (x, ■) - it* (a;) u £ (y,-) - u+(y) 



2p 



<C\x-y\^-^ 2 {\ + ^'). (78) 



the residual process - ■= — satisfies Condition (ii) of Theorem [5] 

V e 

with the exponents ft — 2p and 5 — p/2 — 3/2. 



(79) 



Choosing p such that /3 > and S > (it suffices to choose p > 3), and collecting and (|79l) . we see that the 
random process - -j= satisfies the assumptions of Theorem [SJ It thus converges in law to the Gaussian 



v/i 

process Gq(x,uj) defined by (|2"2")l . This concludes the proof of Theorem^] 



□ 
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5 Technical proofs 

We collect here the proofs of Lemmas [U [TU] and [11] 



5.1 Proof of Lemma 

Lemma [5] is a consequence of the following result: 

Lemma 12. Assume that d<j> satisfy ©, ©, © and ([8]). Assume furthermore the independence conditions |19 

and (|2~lj) . For any < a < (3 < 1, define the random variable 



1 



(80) 



where the function ip is defined by (|17[) . For any p € N*, there exists a deterministic constant C p independent of e, a 
and fj such that 

Ve > 0, E [Z £ (a, [3, -) 2p ] < C p \(J3 - af + e^ 1 )/ 2 " . 
We first prove Lemma [T2l and next Lemma [5J 



Proof of Lemma QH Using the variable s = 



co , we write 



Z s (a,f3,u) 



(it = Vff 



* _1 (/3/e^) 



4> 1 (a/s,Lu) 



ip(s)(j) , (s^ uj)ds 



(81) 



For future use, we introduce, for any x £ (0, 1), the notation 

K x (cj) = 14>-\x/s,u)\. 

In view of dU) and ©, we have 



f3-a 



< 



a 



/3 — a 



Hence, up to some boundary terms (due to the fact that (f> (a/s, uj) and <f>~ (/3/e, to) are not integer numbers), Z e 
is a sum of the variables Yfc defined by (p~5|) . with a number of terms of the order of e _1 . Note however that this number 
of terms, equal to Kp(w) — K a (uj), is random. To proceed, we write Z e as the sum of two contributions: (i) a sum of 
the variables Yfe with a deterministic number of terms, and (ii) a remainder, that will be successively estimated. 
Following (|81D . we have 



ip(t)(j)'(t,U])dt 







ip(t)(f>'(t,u})dt + 



_1 (a/£,cj) 



il)(t)cj)'{t,uj)dt 



= B e (a,/3,u) + A e (f3,uj) - A e {a,ui) 



with 



AJx,w) 



i/}(t)(f>'(t,L))dt, 



■4>(t)(/>'(t,uj)dt. 



(82) 

(83) 
(84) 



Note that, up to boundary terms, B e (a, f3,ui)/y/e is a sum of the variables with a deterministic number of terms. 
We infer from (|82|) that, for any peN*, 



E[Z E (a,(3,-) 2 P] <C P E [B E (a,f3,-fP] + C P E [A e {a, -) 2p ] + C P E [A e ((3, ■ 



|2pl 



(85) 
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where the constant C p only depends on p. We now estimate B £ , and next A £ 
Step 1: Estimation of B £ 



Denoting by K a 



a 



and K 



P 



P 



.eE(/ V). 



K B -X 



fe+1 



.eE(/ V 
B e (a,/3,w) = X] 



, we have 



1+A" 



ip{t)<j)'(t,uj)dt+ r" 1 *' 1 tp(t)4>'(t,u})dt- 



ij}(t)(t>'{t,uj)dt 



k=l+K a 



■>l){t)<j>'{t,u)dt + yfe 



1 + K a 



ip(t)<f)'(t,u))dt, 



(86) 



where we recall that is defined by (TT51) . We thus obtain, for a deterministic constant C p that only depends on p, 

2p 



|B e (a,/3,a/)|*<C p£ * 



Kg-l 



fc=l+K a 



+ ^11^11*^11^11* 



£=°(Rxn)' 



(87) 



Recall that (Yfc)fcez is a sequence of independent identically distributed variables, with E(Yfc) = 0. We now use the 
fact that any such variables satisfy the following bounds: 



E 



V fc=l / 



2,/ 



< 



Cp 
~Np 



for a constant C p that depends on p and the moments of Yj~ , up to order 2p. This is proved by developing the power 
2p of the sum, and then using the fact that the variables are i.i.d and have mean value zero. In our case, the variables 
Yfc are bounded almost surely, and thus all their moments are finite. We thus deduce from (l87l) and ([88]) that 



E[B e (a,f3,-) 2p ] < C p e p E 



Kg-l 



£ * 



k=l+K a 



2p 



+ C p e^\\l% m U'\\ll 



L°°(Rxn) 



< C v er{Kp - *T a - If + C p e?U\\% m U'WT-^.n) 

< Cp {P-aY + C p e^\\% m U'\\% 



l»(lx!l)' 



(89) 



Step 2: Estimation of A £ 

We now bound A e (x, w), for any x e [0, 1]. Denoting K x = 



hence 



thus 



_eE(/ V)_ 



and i^x(^) = w)J, we have 



A £ (x, uj) — \fe | / ip(t)4>'(t,uj)dt+ / i/}(t)<f>'(t,w)dt+ _ ip(t)<f>'(t,w)dt\ , 



|A:(a:,w)| < 



ijj(t)<p'(t,L))dt 



A'.r 



+ 2v / £||V'IU~(r)II0'IIl~(rxo), 



E [A e (x, •)*] < C p e p E 



A'., 



A',, 



4>(t)(f>'(t,-)dt 



2p 



+ C^||V||* (K) ||^||* (RXQ) . 



(90) 
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Let us now bound the first term of the above right-hand side. The difficulty stems from the fact that the random 
variable K x (cj) is not independent from the random process <fi'(t,u). We write, using the bound (|88fl and Young's 

7 

inequality with parameter . 2 +2 > (where 7 > is arbitrary), that 



E 



A'r 



J 



ip(t)<l)'(t, -)dt 



K.r 



2p" 



,/6 



< 



< 



< 



E E 

y- 1 

£—1 OA 



A", 



tp(t)(j)'(t, -)dt 



A", 



2 P 



E 



2j2p+2 



A x +j 



ip(t)cl)'(t, -)dt 



A, 



■2p+2 
' 1 

27 



fc(^) =K x +j] 



7 



rE 



/ 2j2p+2 

7 , J 



E C *P ^2 + 



E n 



fe=A x 
2p+2 



4p" 



■2p+2 



2j 2 2 7 



< c 2p 2+i-E 



|2p+2 



2 7 

K x {u)=K x +j 

I 2p+2 



[tf x (w) = K x +j] 



We are now left with bounding from above E ^ | K x — K x | P+ ^ . To this aim, we first bound from above | K x (w) — K x 
\K x {u)-K x \ 2p+2 < C p (\K x -r\x/e,u)\ 2p+2 K x {u)\ 2p+2 ) 

Recall now that, in view of (O, we have \a — b\ < u~ 1 \cj>{a, u) — <j){b,w)\ for any a and 6, almost surely. We get 



(91) 

2p+2 



\K x {l))-K x \ < 



j/2p+2 
fe+1 



4>(K x ,u) 

£ 



2p+2 



,2p+2 



(92) 



We now recall that the random variables -Dfe(w) = / 0'(t, w)df, introduced in (|2U)) . are assumed to be i.i.d. random 
variables. Writing (f>{x,u)) — 0(0, w) + / (f>'(t,(jj)dt, we obtain that 



(f>(K x ,uj) 

£ 



2p+2 



< a, 



< c n 



A,, 



^'(t,o;)<ft- Jf a E(Do) 



2p+2 



+ |0(O,c)| 2p+2 + M(I? )-- 



2p+2 



A x ~l 



2 (D fc H-E(£> )) 



2p+2 



|0(O )W )| 



2p+2 



2p+2 



where, we recall, = 



Observing that E(D ) = E 0'^ , we have - - A' X E(D ) < E(D ), thus 





2p+2 / 


4>(K x ,uj) - - 


< C P I 





Ka>-1 



2 (-Dfe(w) - E(D )) 



fc=0 



2p+2 



|0(O,o;)| 2j,+2 + |E(£)o)| 2p+2 



(93) 



Collecting and ([55]). we obtain 



|«= / A rr | 2 P+ 2 / ^P 

^(w) - K x \ <-n£ 



A x -1 



£ (D k (u)-E(D Q )) 



2p+2 



+ |<K0, w )| 2p+2 + |E(A>)|^ + z/ 



2p+2 1 ,,2p+2 
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Next, we take the expectation of the above inequality and use ([55)1 to get 



a 



Since K x 



eE(/ V). 



e(\K x - K x \ 2p+2 ) < -gj (i^ + E (WO, -)| 2P+2 ) + |K(Ub)| 



we know that K x is of the order of 1/e, and thus 



2 P+ 2 -|_ j/ 2 P+ 2 



(0,1), e(\K x (u) 



|2p+2 



< c, 



p £ p+l 



(94) 



for a constant C p independent of e and a;. We infer from (j9"Tj) and (l94l) that 

2p' 

V7 > 0, E 



ip(t)<f>'(t,-)dt 



A',,. 



7 



Taking 7- 1 = e^ 1 )/ 2 leads to 



E 



tp{t)(j)'{t,uj)dt 



A,, 



2p" 



< c, 



1 



p £ (p+l)/2- 



Collecting (J5DJ) and ([55]). we obtain 

V.t G (0, 1), E [A s (ai, -) 2p ] ^C^- 1 )/ 2 
for a constant C p independent of e and x. 

Step 3: Conclusion 

Collecting (|55j) . and (l§6l (which is legitimate since 0<a</?<l),we obtain 

E [Z e (a, p, -) 2p ] < C p Uf3 - a) p + e p + e {p - 1)/2 \ < C p \{fi - a) p + £ ^ 1]/2 
where C p is a deterministic constant independent from a, /3 and e. This concludes the proof of Lemma \T2~\ 



(95) 



(96) 



□ 



Proof of Lemma\^ The result directly follows from Lemma [T2l and an integration by part argument. We consider the 
random variable Z £ (a,(3,uj) defined by 



dt. 



Integrating by part, we see that 

Z £ (a,j3,oj) = [A(t)Z £ (a,t,u)f a - / A'(t)Z e (a,t,u)dt = A(P)Z e (a,P,u)- / A'(t)Z e (a,t,u)dt. 

J a J a 

Using the Cauchy- Schwartz inequality, we obtain 

Z e {a,(3,uj) 2 < 2A(P) 2 Z £ (a,/3,uj) 2 + 2 / (A'(t))dt Z £ (a,t,uj) 2 dt. 

J a J a 

We now take the power p of this estimate: 



(a,/3,w) 2p < C p \\A\\f^ atP) Z e (a^M' p + C P \\AT3 {a ,0) [ I Z e {a,t,u}fdt 



,2p 



2p 



,/i|2p 



< ^PII^ (Q ^^(«,/3^) 2p + C P P'||f 2(Qj3) / Z £ (a,t,u) 2p dt 



p/q 



dt 



< 



^MII^^^K/J^^ + Cp^-a^MII 2 ^^ / Z £ (a,t,u) 2p dt 
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where we have used Holder inequality with 1 = 1/p + l/q, Using Lemma IT2l we thus obtain 



E[Z £ (a,/3,.) 2p ] < CpMfc a fl ) E[Z e (a J AO 2p ]+Cp08-a) p -Ml^ / ||? a(a fl ) / ^ [Z £ (a, t, -) 2p ] dt 



< C P \\A\\ 



3 («,/3) 
L~(a,/3) 



{(3 _ a) p + e (P-D/a + _ a) P-i H^'H* / (t - a)" + e^' 2 



r 



< ^PH 2 L~( Q ,« (^-a) p + ^/ 2 +C p (/3-«r 1 p'llf 2(a , /3) (/?-a) p+1 + (/3-a) £ (p - 1)/2 



< C v 



(/3 - a)" + e^' 2 \\A\\% {a ^ + ((3- af \\A'\\f 2(c 



This concludes the proof of Lemma [5] 
5.2 Proof of Lemma [10] 



□ 



By definition, 



Z e (a,/3,u}) 



f A(t)ij) y<t>-' y-,u ) ) dt. 



We start by replacing the function A by a piecewise constant function A, that we will choose later as an accurate 
approximation of A, in a sense to be made precise. We thus introduce the function A defined by 



A' 



A(i) = Y,A\t p ,t p+l) {t), 
p=l 



(97) 



with a = t% < t% < ■ ■ ■ < ijv+i = P- Hence the sets (t p ,f p +i) are disjoint one from another, and Ui< p <at [t p , t p +{\ 
[a, /?]. We associate to this function A the random variable 



(98) 



Step 1: Z s (a,j3,u)) converges in law to a Gaussian random variable 
In view of (|5D1 and (|82p. we have, for each p, 



t 



I i>\4> 1 \j?-> W ) ) dt = Z e(tp,t p+ i,Uj) = B £ (tp,t p+ i, LO) + A £ (t p+ i,Uj) - Ae(t p ,u). 



We can write B £ (see (|86p) as 
where 



B e (t p ,t p+ i,uj) — B £ (t p ,t p+ i,uj) + i? £! p(w), 
B e (tp,t p+1 ,uj) = yfe ^2 Y k{u) 

k=l+K p 



(99) 



with i-C, = 



and 



'p+i 

• / eE(/o *° V(*)0'(t,w)dt- 



_- / ip(t)(f>'(t,w)dt. 



We hence write 



AT 

Z e (a,P,w) = J^-4 P {B e (tp,t p+1 , uj) + R £ , p (uj) + A £ (t p+1 ,u) - A £ (t p ,uj)J . 



(100) 
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Observe that R StP satisfies 

|-Re,pM| < 2v / £||'0l|L~(R)||0'|U'»(Rxn), 

and hence goes to as e — > almost surely. 

In the sequel, we first show that A £ converges in probability and thus in law to as e goes to 0, and next that B e 
converges in law to a Gaussian random variable as e goes to 0. 

Step la: A £ (x,uj) converges in probability to 

For any x € [0, 1], we have 

»(j)~ 1 (x/s,u) 



A e (x,uj) = \fe 



ip(t)c/)'(t,u))dt- 



tp(t)(j)'(t,uj)dt 



>K x {u)-l K x -1 
fc=0 fc=0 



(101) 



where K x 



and = [(j) 1 {x/e,oS)\, and 



l (ar/e,w) 



ip(t)(j>'(t,u})dt 



Jf*(u>) 



ip{t)(j>'(t,u)dt. 



K, 



We have 

(KxQ)j 

hence i? e goes to as e — ► almost surely. 
For any A > and S > 0, we write 



(102) 



K*(w)-1 Ka-l 

E w- e y *n 

fc=0 fc=0 

F„(w)-i a x -i 

E r»( W )- £ nn 

fe=0 fc=0 

A' -1 



> 



> -= and li^xM - -Ka| < 



E w - E w 

fc=0 fc=0 



> -= and |i^xM - -RTjb > 
v£ 



s 


) 


e 





Remark that 



Jf,(w)-i 

E w 

fe=0 



E y *M 



fc=0 



> — p and Iff^u;) — K x \ > 
v £ 



<P Jf x («) - KJ > 



We next write, using that is a sequence of independent identically distributed variables, that 



K~-l 



E y *M - E w 

fc=0 k=0 



> —= and |if x (£j) - K x \ < 
v £ 



< 



sup 

l<fc<[<5/ej 



E y H 



> 



(103) 



(104) 



We now recall the Kolmogorov inequality |HJ p 175]: as Yi is a sequence of i.i.d. random variables with mean zero, we 
have 

k 



sup 

Kk<n 



E y «M 



> x < x _2 Var E y M = na; ~ 2 Var ( r o) • 



We thus deduce that 



K.M-l A x -1 

E w- E y fc( w ) 

fe=0 fc=0 



> -= and \K x (uj) - K x \ < 
v£ 



< ^Var(y ). 



(105) 
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Collecting (f!D5]l . (jTIHf and (fTU5|) . we obtain 



K~-l 



E n( W ) - 2 r ^ 

fe=0 fc=0 



> -= I <p( I^M-^I > 



-^Var(Y ). 



For any fixed A > and any rj, we choose <5 > such that <5Var(Y"o)/A 2 < Recall now that e(f) 1 (x/s,uj) converges 
to j as e — > a.s. (see |7J), which implies that P (|/\z(u;) — if x | > [|J) goes to when e — > 0. There thus 

E (Jo V) _ 
exists £o such that, for any e < Eq, we have P (\K x (lj) — K x \ > [|J ) < 77/2, and thus 

A 



K~-\ 



E w- E y ^ 

fe=0 fe=0 



We thus have proved that, for any A > 0, we have 



lim I 

e-s-0 



if x (w)-i 

E w- E ^-h 

fe=0 fc=0 



> A = 0. 



Collecting this limit with (|101[) and (| 102[) . we obtain that A e (x,oj) converges in probability to as e — > 0, for any x. 

Step lb: Convergence of B e and of Z e 

Recall that (Yk(u>))k£Z is a sequence of i.i.d. variables of mean zero (see assumption (Jl9j) ) - Using the Central Limit 
Theorem, we obtain that B e (t p ,tp+i,u>) defined by (|99p converges in law to a Gaussian variable, 



B s (t p ,tp + i,uj) —^7V(0,fj p ), 



the variance of which is 



Var(y )- 



p mi *o 

In addition, the random variables B e (t p ,t p+ i,uj) are independent one from another. 

As R e (ui) and A s (x, w) converge to zero in probability for any x, and B e (t pi t p+ i, uj) converges in law for any p, we 
deduce from (|100p that Z E (a,/3,uj) converges in law to a Gaussian variable, 



Z e (a, (3, u) Z {a, /3, u) ~ JV(0, ?(a, /?)), 



the variance of which is 



Var(F ) 



JV N 

a{a, P ) =E4^ = E^ iTTT^f Var(r ) - 

p=i E (Jo </> ) E (Jo ^ ) 



.4 



L 2 (a,/3) 



Step 2: Convergence of the random variable Z e (a, j3, u) 

Recall that A is piecewise continuous with a finite number of discontinuities located at {ifc}i<i.< m and that, for 
each 1 < fc < m — 1, A' £ L 1 ^, tfc+i). Introduce the broken Z^-norm of A': 



\A'W{a,(l) ■= MV(a,ti) + E W A 'hHt k ,t k+1 ) + \\A'\\ L i {tm jy 



fc=l 



Let us fix some rj > 0, and let us complement the previous set of points (t p )i<p<Ar+i such that 

a = t\, ijv+i = P and < t p+ i — t p < i] for any p. 



(106) 
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We set 



A p — A(t +1 ) 



and consider the function A and the random variable Z e (a,(3,u>) defined by ([57]) and (|98p. 
We write, for any £ g R, 



JtZ (a,j3,-) 



where Za(a,/3, u>) is a Gaussian random variable distributed according to J\f(0,a(a,/3)), with the variance 

_. „, Var(Fo) „ , ll2 

We successively estimate the three terms of the right-hand side of (|108p . 
For the first term, we first see that 



We next compute 

Z 



< E 



(a,p,w)-Z e (a,p,u) = ±=J (^(t)-^(t)W(W~,w) 



Using the random variable 



Z e (t p ,x,w) = —= ip ( 1 ( -, w ] ) rft, 



we write 



Z e (a,p,u)-Z e (a,P,u>) = (A(t) - A p ) 

* — i ■'tm 



dZ £ (t p , t, lj) 
db 



dt 



^2l(A(t) - A p )Z e (t p ,t,cj)] t l +1 -J2 / A'(t)Z e (t p ,t,uj)dt 
p=i p p=i"'*p 

~zJ / A'(i)Z e (t p ,t,w)dt 
p=i •'h 



where we have used (|107[) . We thus have, using Lemma [T^l that 



E 



Z e {a,P,-)-Z e {a,P,-) < V / |^'(t)| E |Z e (*p, t, -)| * 



p=i " 

w /-Vn ma 

< \A\t)\ ((t-t p y + ^y' i dt 

P =i •' t p 



N 



1/4 



|.4'(t)| dt 



(107) 



+ E / e «W,)_ e ^»W,')) i (i 08 ) 



<|£|E Z s (a,P,-)-Z E {a,P,-) . (109) 
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where C is a constant independent of e and {t p )i< p <N+i- In view of (|106p . we have 



E 



Z e (a, f3, ■) - Z E (a, (3,-) <C (r, 2 + J!) 1 '* \A'\ LH 



Inserting (|110[) in (|109[) . we deduce that, for any e and 77, 



' e ^Z e (a,fi,-) _ e i£Z s ( a ,f3,-) 



(110) 



(111) 



We next turn to the second term of the right-hand side of (|108p . We recall that Z £ (a,(3,uj) and a(a,f3) depend 
on 77, through the choice of the function A. For the parameter 77 that we have chosen, Z e {a, /3, uj) converges in law to 
Zo(a, (3,u)) when e — > 0. Thus, there exists £0(77), that depends on 77 and can be chosen such that £o(f]) < rf'i such 
that, for all e < £o( r ?), 



< 



(112) 



We finally turn to the third term of the right-hand side of (|108D . Since Zo(a, f3,uj) and Zo(a,(3,uj) are Gaussian 
random variables, we see that 

E _ e i(z o (a,0,-)\ = ex p(-^5(a, /3)/2) - exp(-£ 2 a(a, /3)/2). 



Denoting by L the Lipschitz constant of the function a 1— > exp(— £ 2 cr/2) on [0, 00), we thus have 



We next write 



<L\Z(a,P)-a(a,P)\=L^^ 

E(/„ <P') 



A 



~ MIIl 2 (q,/3) 



A 



L 2 (c 



N ftp+i N ft +i 



£ r 

p=i ^ 



In view of (|107|) . we have 

Vt e [t p ,t p+1 ], A(t)=A(t- +1 ) 
Inserting this relation in (|114j) . we obtain 



tp+i 



A'(s)ds = _4 p 



.A'(s) ds. 



A 



Thus, in view of the choice (|106|) . we have 



W A \\h(aJ)=T, f P+ \A+At)) A'(s)d3dt. 



A 



Mllz, 2 (a,/3) 



N ftp+1 ftp+i 

< 2||^|U- (0)W V / / \A'(s)\ dsdt 

p =i ■/* 

A r tp+i 

< 211^11^(^)5^/ |^'(s)| (s-t p )ds 

^ 1 « tn 



w rV; 1 



f p+1 

< 2r,\\A\\ LX{at0) J2 / |^(*)| 

a,/3) • 



(is 



< 



Inserting (|1 15[) in f|113[) . we deduce that 



E(J Q 0') 



(113) 



(114) 



(115) 



(116) 
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Collecting ([TPS]) . (fTTTj) . pT^|) and pTB)) , we have, for any 77 and any e < s (rj) < i] 4 , that 

Var(y ) 



E 



V 7 V ' E(Jo ^ ) 



The above bound holds for any e < £0(77)1 and ?y is arbitrary small. In addition, |^4'|ii( ail g) is independent from 77, even 
though the set of points (tp)i<p</y+i depends on r\. This means that 



lim 



= 0, 



hence Z £ (a, /3,ui) converges in law to Zo(a, /3,w). This concludes the proof of Lemma 1101 

5.3 Proof of Lemma 1111 

According to the definition f|63[) . we have: 

1 F 1 ft \\ , , 1 



\H e (x,u) - H e (y,u)\ < 



dt 



+ \x-y\ 



Thus, for any p £ N* , using Lemma [5] and the fact that \y — x\ < 1, we have 

1 



,oj) )dt 



E 



\H e (x,-)-H e (y, 



,2p 



< C P E 



V? Jy 

+C p \x - y\ 2p E 







2p- 













y*»( f, ( 





2p" 







< C p (\x-y\P + e^/ 2 ) +C p \x-y\ 

< C p (\x-y\P + e^l 2 ) . 



2p 



This concludes the proof of . 

Assume now that \x — y\ < e. We infer from (|117[) that 



\H e {x,uj)-H £ {y,uj)\ < 



x-y 



MU°°(R) (||"4l||i~(0,l) + \\M\\l^(o.i)) < Cy/\x - y\, 



(117) 



where C is a deterministic constant independent of e, x and y. This concludes the proof of (JBSJ), and hence the proof 
of Lemma [TT] 



6 Approximation of the homogenized matrix 

The aim of this section is to prove our second main result, Theorem 0] Since the approach described in Section 13.21 
mimicks the approach proposed in [10] , our proof essentially follows the arguments used in [10] . Because our proof is 
involved, we feel that it is useful to first recall the arguments of [10] in Section 16.11 We then collect some technical 
results in Section 16.21 before turning to the actual proof of Theorem 2] in Sections 16.31 and | 



6.1 Convergence proof in the classical random homogenization setting 

Consider the classical random homogenization problem 

— div 



A[^,L)) Vu £ (x, uj) = f{x) in V, u £ (-, uj) = on dV, 



where V is a bounded open set of R d , f £ L 2 (T>), and A is a stationary matrix in the sense of (|3]), satisfying classical 
coercivity and boundedness properties. The associated homogenized problem is where the homogenized matrix is 
given by 



A (y, •) (e 3 + Vw ej (y, •)) dy 
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where Q = (0, l) d and where, for all p £ M. d , w p solves the corrector problem 

-div [A(y, u) (p + Vw p (y, w))] = in R d , 

Vw p is stationary in the sense of ([3]), E (^J Vw p {y, )dy^j = 0. 

In [10) . the following approximation strategy is proposed: introduce the approximate corrector w p (-,u>) as the Qn~ 
periodic function satisfying: 

for all ^ Qw-periodic, f (Vi/>) T A{-, oj) (p + Vwg(; u)) = with f w£(;w)=0 (118) 
and the approximate homogenized matrix A^(uj) defined by, for any 1 < i,j < d, 

[A * N]ij(U})= \ch\I ^(-. w )( e i + V<0,w)). (119) 
Then (see [TOJ Theorem 1]), we have that 

lim A* n {uj) — A* almost surely. (120) 

JV— >oo 

A key ingredient of the proof is the following classical homogenization result (see [T5J Theorem 5.2 p. 151]): 

Theorem 13. Let A E be a sequence of matrices that G-converges to A* in a domain V , and let V\ be an arbitrary 
subdomain ofV. Let p £ W 1 , and assume that the functions w p £ ff 1 (Vi) satisfy the conditions 

w P «£° weakly in H l {Vx), and - div [A £ (p + V^)] =0 in T>'{Vi). 

Then we have that 

A E {p + Vw e p ) A*(p + Vw p °°) weakly in (L 2 (V 1 )) d , 

where u>£° satisfies 

-div[A*(p + Vw™)] = inV'{V x ). 
The proof of (|120p goes as follows (see [10] for details) . The rescaled corrector 

Wo )P (x,uj) := -^w* (Nx,u) 

is shown to satisfy the a priori bound \\wq p (-, uj)\\h 1 (q) < C, where C is a deterministic constant independent from N. 
We thus deduce that, almost surely, there exists a Q-periodic function w^ p (-,uj) £ H 1 (Q) such that 

Wo!p(->w) w™ p (;u) weakly in H l {Q). 
Consider a Q-periodic function i/j £ H 1 (Q). Choosing "4>n{v) = ^(y/N) as test function in (I118p . we obtain 

(VlPf A(N;L0)(p+Vw£ p (;Uj)) =0. (121) 

We are then in position to use Theorem fl3l on the domain V\ = Q. We thus get that A(N-,uj)(p + Vw^) weakly 
converges to A*(p + Vuig° p ) in (L 2 (Q)) d . We then infer from (|f 21D that, for any Q-periodic function ijj, we have 



(V^) A*( P + Vw% p (-,u)) =0. (122) 
This implies that Vw^° p (-,w) = 0. Using the same weak L 2 convergence as above, we deduce from (|119p that 

[^]« H = j Q ejA{N-,u) (e 3 -+V< ej (., W )) ^efA* (e, + V< e , (•,*)) = [A% . 
This concludes the proof of (|120l) . 
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6.2 Some technical ingredients for our analysis 

A key ingredient to prove Theorem 0] is to find an appropriate domain on which to apply Theorem 1131 The following 
lemmas are useful for that purpose. 

We first recall (see [7J Lemme 2.1]) that — <fi(N-,uj) converges to E (^J V</j ■ in Lj^ c (M d ) almost surely. Likewise, 



Q 

in view of the proof of [7J Lemme 2.2], we have that — <j)~' 1 {N-,ui) converges to 
surely. The functions being smooth, we thus have that, for any compact K, 



E ( / V0 

Q 



in L^ c (R d ) almost 



lim 

TV— too 



l#JV-,w)-E| / Y<> 



= lim 

iV-KX 



l0-i(JV-,w)- 



E / Wcj) 



C°(K) 

As pointed out in the proof of [7J Lemme 2.2], a consequence of the above fact is that 

= a.s. 



= a.s. 



C°(K) 



lim 

TV— >oo 



(123) 



(124) 



This can be shown by first assuming that </>(0, ui) — 0, and using a regularization of the indicator functions. The general 
case <fc(0,ui) ^ next follows as an easy consequence. 

The first ingredient we need to prove Theorem 2] is the following lemma, which is somewhat related with the above 
results: 

Lemma 14. Let <j> be a diffeomorphism that satisfies (J7J) and ([8]). For any compact set K that is a proper subset 
of the open set E (^J V</>^ Q, and almost all uj, there exists Nq[lo) £ N such that 

viv>TV M, k c^(Qtv,w), 

o 

where K denotes the interior of the set K and, we recall, Qtv — N Q. 

The following easy result is useful to prove Lemma H4l 

Lemma 15. Let <f> be a diffeomorphism that satisfies ([5]) and (J7J. Then there exists a deterministic constant LLip such 
that the diffeomorphism _1 (-,u;) is Lipschitz with that constant. 

Proof of Lemma\15\ We infer from (J7J that V(/; T V0, which is a symmetric matrix and therefore diagonalizable, has a 
bounded spectrum. The assumption ^ then implies that the eigenvalues of V0 T V</> are bounded away from 0. Hence, 
there exists a deterministic constant c > such that for all £ € R d we have 

£ T (V(t)(x,uj) T V(f>(x,Lj)£, > c\£\ 2 a.s., a.e. on R d . 

For any £ G R d , we set £ = (V0(x, uj))^ 1 ^ and obtain that 

\(VcP(x,cj))-^\ Kc- 1 ' 2 (125) 

The diffeomorphism iT; -1 (•, oj) is thus Lipshitz with the deterministic constant c~ x / 2 . □ 

Proof of Lemma^TJ^ Let if be a proper subset of the open set E (^J V</3^ Q, and let us fix uj such that 

1 



— 4> (N-, uj) converges to 



E / V<f> 



in C°(K). 



(126) 



In view of p23[) . we know that (1126[) holds for almost all uj. 
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We prove Lemma [TJ] by contradiction. Suppose that, for all iVo G N, there exists N(Nq,lj) > Nq such that K is 
not included in „, „ r<KQiv(iVo,w)> w )' Otherwise stated, there exist N(Nq,uj) and z(Nq,uj) such that 



N(N a ,u) 



z(No,u) € K and z(N ,u) 



1 



N(N ,lo) 



<l>{N{No,u)Q,u) 



Introduce y(N ,uj) = 



1 



AT(iV ,a;) 



4>- 1 (N(N ,uj)z(No 1 uj),uj). We thus have that 
y(N ,oj)^Q. 



(127) 



We now pass to the limit No —> oo. Observing that z(No,u>) belongs to the compact set K, we deduce that 
{z(Nq, w)} n gN is a bounded sequence and thus converges, up to the extraction of a subsequence, toward some ~z{uS) G K. 

Let us now show that {y(N ,uj)} No<£N is also a bounded sequence. Using the fact that the diffeomorphism w) 
is a Lipschitz mapping with a deterministic constant Lup (see Lemma fT5|) . we write 



l»(JVo,w)| 



N(N a ,u) 



U- x (N(N Q ,u>)z(N ,u)),u>)\ < L Li M N o,u)\ 



N(N ,u) 



We deduce that, almost surely, {y(No,uj)} No£N is a bounded sequence and thus converges, up to the extraction of a 
subsequence, toward some y(co). In view of (I127[) . and since Q is an open set, we have that y(u>) ^ Q. 
We now claim that 



(128) 



Indeed, we write that 











v(N ,u)- 






z(w) 



< 



< 



1 



N(N ,u) 
1 



0- x (iV(iVo,u;>(iVo,a;),a;) 



E( / W(f> 
Q 



z(N ,lj) 



z(N ,u) 



E( / V</> 
Q 



N(N ,u) 



+ C?|*(JVo,w) -«(«)| 



C°(if) 



Both terms converge to when No — > oo, respectively in view of (|126D and of the definition of z(lu). By definition of 
y {<*>), we deduce (|128l) . 

We now reach a contradiction since ^(w) E K C K Wcj^j Q whereas y(uj) £ Q. This concludes the proof of 
Lemma [TU 



□ 



The second ingredient we need to prove Theorem [4] is the following lemma: 



Lemma 16. Let <f> be a diffeomorphism that satisfies ^ and ([7]). There exists an open set Q{uj) and some k{u) G N 
such that 



E / V(j) QcQ(cj), 



ViV g N*, QH c -<b(QHu)N,u). 



(129) 
(130) 
(131) 



Proof. The first assertion relies on the fact that, in view of (0, we have 

ViV G N*, —\(f>(Nx,cj)\ < M\x\ + -^^((Lw)! < M +\<j)- 1 (0,Lj)\ a.s., a.e. on Q. 
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It is thus sufficient to choose Q(u)) such that 



— M — 



L (0,c<j)|,Af + \(f> 1 (0,w)| C Q(u}). Upon choosing a larger 



Q(lo), the second assertion is also satisfied. Now that Q(to) is chosen, we show that we can choose fc(w) such that the 
third assertion is satisfied. Using Lemma [15l we see that, almost surely, 



1 



1 



ViVeN*, -M-^Nx^l <Lu P \x\ + -\<K0,w)\ < L Up N + |^(0,w)| a.e 



N 



N 



There thus exists k(u) such that, for any N £ N*, we have —(f) 1 (jSfQ(u), u?j c Qk(u)- This implies the third assertion 



and concludes the proof. 

6.3 Proof of Theorem H 

To simplify the notation, we introduce the matrix 



□ 



E 



Q 



adxd 



As pointed out in [SJ Remark 1.9], we have that 

dct a = E 

We hence deduce from © that 



det 



(132) 



(133) 



det a > v > 0. 



We also introduce the matrix (3 g 



pdxd 



defined by 



P = E 



HQ,- 



= E 









det(V(/)) (V0)" 



(134) 



The proof of the following lemma, useful for proving Theorem 21 is postponed until Section 

Lemma 17. The constant matrix (3A*a~ T is coercive. 

The proof of Theorem|3]is composed of four steps. In Step 1, we introduce a rescaled corrector, denoted Wq p (-,v) 
(see (|136[) below), and show that it converges toward some function uig° p (-, lu) weakly in H 1 . Then, in Step 2, we prove 
that w™ p (-,w) is aQ-periodic. Next, in Step 3, we show that w^ p (-,oj) solves the equation — div [BVw^,] = in R d 
for a constant deterministic matrix B (see (|151[) below for a precise statement). Combining these results and using 
Lemma IT71 we conclude that Vw§° p = 0. This is a key ingredient to prove, in Step 4, that the random approximation 
A* n (uj) indeed converges to the homogenized matrix A* almost surely. 

Step 1: Introduction of a rescaled corrector Wq p , and convergence of Wq p to some uig° p 

We first establish some a priori bounds. Taking ijj = w p as test function in (j2"4"]l . and using ^ and |7]), we see that 



||(V^(-, w ))- t V<(-,c)|| l2(Qn) < cVlQ^l 
where C is a deterministic constant independent from N. Using again ([7]), we deduce that 

\\Vw»M\\ l2{Qn) <C^/\Q^\. 
Let k € N. Since w p is QAr-periodic, we infer from the above bound that 



\\Vw?M\\» (Qhs) <C,/\Q*n\ 

where C is a deterministic constant independent from N and k. 
We now introduce the following rescaled corrector: 

Wq, p (x,oj) = jjW p (Nx,uj), 



(135) 



(136) 



29 



where, we recall w p (y,u)) = Wp (cj> (y,w),Lj). Using ([7]) and (|125p . we infer from (|135[) that 

where C is a deterministic constant independent from N and We now choose k in the above bound equal to the 
integer k(w) defined in Lemma UHl We infer from the above bound and (|131[) that 



VNeN*, ||V< p || L2((5M) <C(w). 



(137) 



Recall that the solution u;^ to is unique up to an additive constant. We now fix this constant by choosing 

= 0, where the set Q(w) is defined in Lemma 1161 In view of (|136[) . this means that 



w p such that . 

lNQ(u>* 



I Wq p (-,u}) — 0. Using (|137l) and the Poincare-Wirtinger inequality, we deduce that there exists C(w) such that 

ViVeN*, ||< p (-,a;)|| ffl(5(a))) < C(w). 
This implies that, almost surely, there exists wg^, (•,&;) G -ff 1 ^Q(u;)^ such that 

- «>{£,(•, w) weakly in tf 1 (q(w)), (138) 



and, using the Rellich Theorem, that 



w^ p (-,a;) w^ p (;0j) strongly in L (Q{uj) 



(139) 



Step 2: iug° p is aQ-periodic 

We infer' from (fTBTjj) that 



w ( <j>{Ny,u) 
\p [ N ■ 



<A = ±w»{<l>(Ny,u,),u,) = ±w»(Ny,u>). 



Since the function is Qw-periodic, we see that the function y i— > w 



n ( H N v^) 



\ N 



ui I is Q-periodic. Hence, for any 



k G Z d , we have, almost surely, 

K p (ay, w) - < p (a(y + fc), w)] 2 dy < 2 / \w% p (ay, u) - < p ( ^ N ^ U \ 



Q 



dy 



dy. (140) 



We now show that both terms in the above right-hand side converge to when N oo. It is sufficient to consider the 
first term. Let us fix rj > 0. 

We observe that the first term in the above right-hand side satisfies 



oo i \ n ( 4>(Ny,u) \1 2 
w , P {ay,uJ) -w^pl — ,uj i 



dy<2[C Ar H + C 1 Ar H] 



(141) 



where 



w^p{ay,uj)- w QtP I — , u J 



dy, 



Q 



0,p 



<l>(Ny,u) 
N 



n ( <P( N y,u) 



J 0,p 



dy. 
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To show that Cq (cj) converges to 0, we use the fact that the function 



<l>(Ny,u) 
N 



converges to the function ay in 



L^ c (R d ) almost surely (see [7J Lemme 2.1]), and a regularization argument. Since w^ p (-,uj) 6 H 1 ^Q(uj)j , there exists 



such that 



We then write that 
where 



Qn» 

We infer from (lTl2"j) and ([1101) that 

ViVeN*, C7&(w) = 
Likewise, we infer from (|T42l) . (|T29l) and © that 



IK 60 (' ) w)-«C(-.w)ll £ »(fl( U) )<'?- 

w™ P (ay, w ) - (ay, w)] 2 dy, 



(142) 
(143) 



<° (ay, w) - 
</>{Ny,u) 



N 



-,u> 



dy, 



N 



7 °>n N 



Vdet ( 



koTp-^IU 2 (aQ) 



< 



Vdet 



=77. 



V7V € N*, Cfi( W ) < -5=||< " "Cll^^)) < 4* 



N 



(144) 
(145) 

converges to the 



We now turn to Cq[(u}). Using the fact that w^(-,lo) € C°° ^Q(w)^ and that the function 

function ay in L^ c (M. d ) almost surely, we obtain that Cq[(oj) converges to zero as ./V goes to infinity, almost surely. We 
thus can choose N(rj, w) £ N such that 

W>N( V ,lo), C$(w)<rj. 



Collecting (fI35]l . (jTHl) . (11151 and ([135]) . we conclude that 

C^(w) ^ as TV goes to infinity, almost surely. 
We next turn to Cf (w), which is non- negative by definition, and satisfies, using (|129[) and (|139p . 



Q(u 



(146) 
(147) 

as AT — >• oo. 

(148) 



1 ^ (Q «'" > det(V^) K^^"^^^] dy - -H w o°p(-' a; ) - w o V p (-^)l!i 2 (Q ( ^ ) ) 
Collecting (|140l) . (|141l) . (|147l) and (11481) . we deduce that, almost surely, 

VfceZ d , / [u^ p (ay, w) - w^ p (a(y + k), ui)] 2 dy = 0. 

JQ 

The function iu§° p is thus aQ-periodic. 

Step 3: w^° p solves div [.BVwjj^,] = in 2?'(R d ) where B is a constant deterministic matrix 

In the two above steps, we closely followed the proof strategy of [TU] recalled in Section loTTl This Step 3 follows 
a slightly different pattern, and is more involved than the corresponding argument in |10j . which consists in showing 
the weak formulation (I122p . As pointed out above, the difficulty comes from identifying an appropriate domain, 
independent oi N, on which to apply Theorem 1131 To circumvent this difficulty, we work on the entire space M. d , with 
test functions of compact support. 
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Introduce a test function ip £ 2?(R d ), and define the Qjv-periodic function 

My) ■= ^ ( Jj y ~ k 



fcez d 



We note that, for any y € Qn, only a finite number of terms in the above sum do not vanish, and that this number of 
terms only depends on the support of ifr and thus is independent of N . 
Choosing i/jn as test function in ((24]) . we write 

J det(V0(»,o;)) Y, V ^(^f-*) (V0(y^))" 1 A pcr (y)(p+(V0)- T (y, W )V^( 2 /,^)) dy = 0. 
After the change of variable z — 0(y,c<j), we obtain 

that we recast, using the definition (I136|) of u>^ p , as 

VW e N* , ^7^(^=0 a.s, (149) 

fc£Z d 

where 



Vfc e Z d , lim (w) = J£°(u>) a.s., (150) 



We claim that 

Vfc e Z d , 

N-too 



Ij?(u):= ( (Vip (a- 1 z-k)f (3 A* (p + Vw^ p (z,uj)) (det a)- 1 dz, 



where 



where the constant matrices a and (3 are defined by (|132p and (|134[) . 

Assume momentarily that (|150p indeed holds. Then, as the sum in (| 149[> has a finite number of terms, independently 



of N, we can pass to the limit N — > oo and obtain that 

EW=0 a.s. 



which also reads 



Y I {^^{a- l z~k)) T (3A* (p + \Jw^ p {z,uo)){&etay 1 dz = Q. 



Using the aQ-periodicity of the function (shown in the above Step 2), we deduce that 



hp 



(Vtp(z)) 1 f3A*(p + Vw^ p (az,Lo)) dz = 

for all test functions ip G 2?(K rf ). We indeed have shown that 

-dw[(3A*Vw™ p (a-,uj)] =0in2?'(R d ). (151) 

To conclude this Step, we are hence left with showing (|150[) . Formally, this comes from the strong L 1 (R d ) convergence 
of the indicator function 1 1 ^q n u -j towards 1 q q and from the div-curl lemma. We indeed observe that the integrand 
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in/f (w) is the product of (V0(0" 1 (A^z,w),w)) T V^> f ^0~ x (iVz,u;) - fcj with Ape^- 1 ^,")) (p + Vw N p (z,w)). 

We will show in the sequel that the first factor is curl-free, whereas the second factor is divergence free. Using the 
div-curl lemma, this product converges (at least in the sense of distributions) towards the product of the weak limits 
of the two factors, which can be identified. One difficulty to make this argument rigorous is to find a fixed domain 
(independent of N) on which to apply the div-curl lemma. For that purpose, Lemma [Ml is useful. 

We now proceed in details. Let rj > 0, and let O v C O v be two deterministic open sets such that O v is proper 
subset of aQ, O v is a proper subset of O n , and 



aQ\O v 



< V, 



o v \o r 



We then decompose 4 (ui) and I%°(u>) as follows: using (|129[) and (|130[) . we write 



1% (w) = /£» + O), 4» = 4°» + 7^», 



(152) 



(153) 



with 

4"» 
4°°» 



V^(j^cl>- 1 (Nz,cj)-k\\ (V^^- 1 ^, «),«)) 1 A pCT (r 1 (^,w))(p + V< p (z, W )) dz, 
(V?/> (a _1 z - fc)) T /3A* (p+V^fz^)) (deta) -1 dz, 

x (V^^-^iVz,^),^)) -1 A pcr (^-^JVz.w)) (p + V< p ( Z ,w)) dz, 



(W(a _1 z- fc)) T /3A*(p + Vwg° p (z,u;)) (deta) _1 dz. 



To use the div-curl lemma, we need to further decompose l£ (u) and (u>). Introducing a smooth truncation function 



(eD (^) sucn that < £(x) < 1 a.e. and £ = 1 on I( , we write that 



4 V » = 4I» - q v («), 4°r» = / fe » ( W ) - (7~H, 



(154) 



where 



4 W » 



4°f» 



£(z) (V^ (a^z -k)) T f3 A* (p + Vw™(z, w)) (det a)" 1 dz, 



£(z) (V^a -1 *- fc)) T /3A* (p + Vio^ p (z,w)) (deta)" 1 ^. 



We first bound from above C^(w), C~(w), and K^ v {ui). As |f| < 1, we see that 



1/2 



\C^{u)\ < ||VV|U-||V0- 1 (.,w)|| L -,||A por || i oc O v \O v \\p + V< >P (;u)\\ L *(o v \o v ) 
< Cy/¥j \\p + Vw£ p (-,u))\\ l2 (q {u) } 



where, in the second line, we have used (|125l) . (| 152[) and the fact that O v \ O v C O v C aQ C Q(oS) (see (|130p 'l. Now 
using (I137p . we deduce that there exists C{w), independent of t] and N, such that 



ViVeN*, \C? (u)\ < C(u)y/fj. 



(155) 
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Wc likewise obtain that 



\C?(u)\ < C(u>)y/rj and \K% v (u)\ < C{u)y/rj. 



(156) 



Now turning to TZ^^u), we obtain, using similar arguments, that 



K tV {w)\<c 



Using (|137p . a triangle inequality and (|152l) . we deduce that 



i 2 (QM) 



Recall now that, in view of (| 124[) . we have 



lim 

N->oo 



L 2 (QM) 



a.s. 



(157) 



(158) 



We eventually estimate If? (u)) — If** (u>) using the div-curl lemma. The compact O n being a proper subset of aQ, 



we infer from Lemma IT41 that there exists Nq(uj) such that for all N > Nq(u}), O v C —</)(Qn,<jj). We then deduce 
from ([M]) that, for any N > N (u), 

-div[A peI (<t>- 1 (Nz,Lj))(p + Vwg p (z,L>))]=Q mV'(d n y (159) 

Using (1138[) . we can thus apply Theorem IT3l on the domain O v , and obtain that 

Vr {4>-\Nz,u)) (p + Vtu&fow)) -± A* (p + Vw^ p (z,u)) weakly in L 2 (£„) . (160) 

From the proof of [7] Lemme 2.2], we know that — cjT l (N z,oS) strongly converges in L^ c (IR d ) toward oT x z. As, by 
definition, V-0 € C°°(R d ), we obtain that 

(161) 



[jj _1 (iVz, w) - fcj V-0 (a -1 ^ - k) strongly in L" c (M d ). 

Since (V^)^ 1 is stationary, we infer from Lemme 2.2] that 

(V<K0 _1 (iV*,a;),a;)) -1 (deta)" 1 /? weakly-* in X 00 (R d ), (162) 
where the matrix /3 is defined by (|134j) . As O n is a bounded open set of M. d , we deduce from (|16ip and (|162l) that 



^VV' (^^(Nz,^) - k\j (V0(^) _1 (^> w )) w )) 1 ^ (deta)" 1 (Vi/> (q;- 1 ^ - fc)) T /3 weakly in L 2 . 



We eventually note that 



(V^Or 1 ^,^),^)) T V^^0 _1 (iVz,cj)-fc^ is curl-free, 



(163) 
(164) 



as this vector is the gradient of ip 1 (Nz, to) — kj . Collecting (|159l) . (|160[) . (j 163[) and (|164l) . we are in position to 

apply the div-curl lemma (see for instance [T31 Lemma 1.1 p. 4]). We thus obtain that 

V77, lim 7^H=J^M a.s. (165) 

N — — - : ' ' ' 



Collecting (Tl53"l) . (|T54"1) . ([TSS]) . ([TMl) . (|T571) . ([1551) and (IT^5|) . we deduce the claim ([I50]) . This concludes this Step. 
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Step 4: Conclusion 

Collecting the conclusion of Step 2 and (|151[) , we have shown that the function w^ p (-,uj) solves the problem 

-div [/9A*Vto~ p (a-,w)] = in V'(R d ), w^ p (a-,u) is Q-periodic. 

The function g(x,u) = w$f p (ax, to) is thus Q-periodic and satisfies 

-div [p A* a - T Vg{-,uj)] = in V'(R d ). 

We know from Lemma [T71 that the matrix f3 A*a~ T is coercive. The above equation has thus a unique solution (up to 
the addition of a random constant), hence V<? = 0, which implies that Vwg° p = 0. We thus deduce from (|160[) that 



A nor (4>- l (N;u)) (p + Vw£ p (-,u)) A*p as N -> oo, weakly in 1? (Or,). 



1 pcr 



(166) 



We are now in position to prove the convergence of the approximation described in Section 13.21 We infer from (|26p 
that 

[^]« 3 -M=0) + 0), (167) 

with 

= /_ tatfgW*) - W*)) ef A pei {(l>- X (Nx,u)) U + V< ej (*,«)) da:, 

JQ(u)) 



JO, v J 

where we have used that O v C aQ, (IT2"9"ji and ([T50"]) , We deduce from (flgo) that 



hence, in view of (I152p . 



Vry, lim R^(w) = \O v \ [A% . a.s. 

N—>-og J 



lim lim = |«Q| [A*] a.s. 



Turning to i?f ?? (aj), we deduce from (113T[) and (|152p that 



|i?^H|<cM 

hence, in view of (I158[) . 



^(Qk.w) - 1 o, 1 



£ 2 (<3M) 



< C(u) 



L 2 (Q(u)) 



Vv 



lim lim R^ v (uj) — a.s. 



rj-yO N^oo 



Collecting (flFfl) . ((1681 and (fT69)) . we obtain 

lim [B^. («) = |aQ| A* a.s. 

We then deduce from (p?5|) the claimed convergence. This concludes the proof of Theorem [5] 



(168) 



(169) 



6.4 Proof of Lemma 1171 

We first show that 



the homogenized matrix A* defined by (ITU1) is coercive. 



For any p £ M. d , we indeed have 



p T A*p = (deta)^ 1 E ( / 



P T A pci (4)- 1 (y, ■)) (p + Vw p (y, ■)) dy 



(deta)- 1 !} f ( P + Vw p (y,-)f A pcv (<j,~ 1 (y,-))(p + Vw p (y,-))dy 
\J<HQ,-) i 



(170) 
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where the last line is obtained using the arguments presented in the existence proof of Theoreme 1.2]. The matrix 
A peT being coercive (see ©), we deduce that there exists C > such that, for any p £ M. d , 



p T A*p > a_(deta)" 1 E ( / 



C 



E 



<MQ,-) 
p + Vw p {y,-) 



(p + Ww p (y,-)) (p + Ww p (y,-))dy 

2 



(Cauchy-Schwarz inequality) 



> C p T p (using third line of (ITTIO. 

This proves (fTTO)) . 

We now claim that the matrix j3 defined by (|134[) satisfies 

f3 = dot a aT 1 . 



(171) 



This is obvious in dimension d = 1, and also in dimension d — 2, using the explicit formula of the inverse of a 



2x2 matrix. In dimension d > 3, we observe that j3 = E 



adj V0 



where adj V<^ is the adjugate matrix (i.e. 



the transpose of the matrix of cofactors) of V0. The matrix V</> being stationary, we deduce from [TJJ Corollary 1] 



adj Vcj) 



= adj 



V<f> 



, from which 



and (|133j) (see also [U Corollary 6.2.2] for the specific case d = 3) that E 
we readily infer (|171l) . 

We are now in position to prove Lemma |T7j Using (|171[) and (|170p . we indeed see that there exists C > such 
that, for any p £R d , we have 

p T (3A*a~ T p = deta p T a~ 1 A*a~ T p > Cp T a~ 1 a~ T p. 
Since deta > 0, we see that the matrix a _1 a _T is symmetrix positive definite, which concludes the proof of Lemma 1 171 
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